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Abstract 

A general self-consistency approach allows a thorough treatment of the correc- 
tions to the standard mean- field approximation (MFA). The natural extension of 
standard MFA with the help of a cumulant expansion leads to a new point of view 
'sj- ' on the effective field theories. The proposed approach can be used for a systematic 

. treatment of fiuctuation effects of various length scales and, perhaps, for the devel- 

I opment of a new coarse graining procedure. We outline and justify our method by 

some preliminary calculations. Concrete results are given for the critical tempera- 
D . ture and the Landau parameters of the (/>^-theory - the field counterpart of the Ising 

Q I model. An important unresolved problem of the modern theory of phase transitions 

- the problem for the calculation of the true critical temperature, is considered 
within the framework of the present approach. A comprehensive description of the 
ground state properties of many-body systems is also demonstrated. 
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O ' 1. Introduction 



. This investigation is focused on the correspondence between microscopic models of phase 
transitions and their quasi-macroscopic (field-theoretic) counterparts. Here we shall out- 
Qj [ line a new self-consistency approach to a more accurate derivation of effective field theories 
' from microscopic models defined on lattices. 



Our method is general and can be used for a wide class of microscopic models but for a 
concreteness here we shall illustrate our approach with the Ising model (IM), given by 

1 ^ 

= --^JijSiSj , (1) 

where s = {si} denotes a lattice "field", Sj = ±1, and the interaction constant of ferro- 
magnetic type Jij = J{\i — j\) > depends on the intersite distance |i — j| in a regular 
/^-dimensional lattice of A^ sites ("spins" or pseudo-spins). Note, that J(0) = Ju = 0. 

We shall follow the main path of the phase transition theory, where the effective (quasi- 
macroscopic) field Hamiltonians (alias Ginzburg-Landau (GL) free energies) are derived 
with the help of two systematic methods: (i) Hubbard-Stratonovich transformations 
(HST) and, (ii) a mean- field (MF) like procedure IH El Ej- Here we propose a more 
thorough approach based on a convenient generalization of (ii). 

^Email: uzun@issp.bas.bg. This work has been partly written during a visit in the Abdus Salam 
International Centre for Theoretical Physics, Triest. 
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The known field theories exhibit both success and failure in the description of many-body 
systems defined by microscopic models. For example, we believe that the renormalization 
group methods of the modern theory of phase transitions PP based on the 0^-theory yield a 
quite convenient description of the scaling and universality properties of IM but we cannot 
be certain that the field theory satisfactory describes important non-universal properties, 
such as, for example, the critical temperature Tc and the lower critical dimensionality Di 
(for IM, Dl = 1, whereas within the 0^-theory, = 2). In fact the field theory fails 
along this line of studies. For example, the fluctuation shift ( AT^) f of Tc predicted within 
the one-loop approximation for the 0^-theory is a very small and, hence, unrealistic, while 
in the higher orders of the loop expansion this shift turns out infinite and its calculation 
needs a special renormalization. As a result, the problem for the value of within the 
framework of the present field theory of phase transitions remains unresolved. Here we 
shall show the genesis of this problem and present a satisfactory solution. 

The reason for the mentioned difficulties of the field theory approach is in the quite sim- 
plified coarse-graining procedures used in the derivation of the GL effective Hamiltonians 
from microscopic models. HST is applied together with the long-wavelength approxima- 
tion (LWLA), (kao) <ti n; k = \k\; k = {ki,...,k£)) is the wave number, Oq is either the 
lattice constant or, generally, the mean interparticle distance. LWLA leads to a correct 
expression for the Ornstein-Zernicke correlation function but cuts the short-range (i.e. 
high-energy: e{k) ~ fc^,A < < 0/ao) interparticle correlations of fluctuation type that 
have the main contribution to the shift (ATc) /; A <C vr/ao is the upper cutoff for k within 
LWLA. 

In the approach (ii) we usually say that we neglect the fluctuations of the physical quan- 
tities from their equilibrium values. This is not entirely true. It seems important to 
emphasize that the latter are values in MF approximation (MFA) and, hence, they are in- 
correct. Thus the fluctuations are deflned towards incorrect statistical averages, and their 
contribution to the free energy of the system cannot be accepted as an entire fluctuation 
effect. The latter can be correctly evaluated, if we are able to deflne the fluctuations as 
variations towards exact statistical averages calculated by the Hamiltonian (^. This task 
seems unsolvable, but the present paper makes a step of improvement of the theory along 
the same direction. 

Note, that the methods (i) and (ii), no matter of the difference between them, lead 
to the same GL effective fleld theory. Both methods use LWLA. For the method (ii) 
LWLA seems to me obligatory because of the following important argument. HST can 
be applied only to positively deflnite matrices (Jij) but this condition is not satisfled 
by quite important lattice models, such as the nearest-neighbour (nn) IM. In LWLA, 
however, the same interaction matrix Jij-nn of IM is modified to a form that is positively 
definite. Besides, in both methods, LWLA is used to help the derivation of the quasi- 
macroscopic (fiuctuation) and macroscopic (thermodynamic) properties, i.e., the LWLA 
is a tool of a "coarse-graining" procedure for the many-body system. But in both cases (i) 
and (ii), the micro- and mesoscopic interparticle correlations are ignored. Here we shall 
use a generalization of the method (ii) in order to improve this disadvantage of the theory. 
In this way we shall present more accurate calculation of the LG parameters (vertices) 
of the effective field theory of Ising systems. Besides we shall show explicitly for the 
first time the mechanism of "statistical correlation" which leads from the two-site {i — j) 
trivial correlation presented by the initial interaction Jjj to short-, meso-, and large-scale 
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effective multiparticle correlations of fluctuation type - fluctuation correlations. Thus 
we shall establish and develop for the first time a new "coarse-graining" procedure, and 
this is the main aim of our report. Besides, we shall show that the GL parameters of 
the 0^-theory acquire (l/2r)-corrections for both short-range and long-range interactions 
J{\i — j\) pp. These corrections will be presented to second order in (1/z)^. 

The 1 / z-expansion has been introduced by R. Brout (4ij and applied in calculations of Tc 
and thermodynamic susceptibilities El 13 IHl El HHl HH 1121 UB] ■ We are not aware of 
another relevant work along this line of research except for our recent investigation [21 E] ; 
the latter will be used in our investigation. Remember, that the (Brout) approach is a 
development of an older method - the Kirkwood method of cumulants (semi-invariants); 
hence, we shall follow a cumulant expansion which is well known. 

We shall generalize the Brout approach in a way that makes possible to derive effective 
field theories, and this will allow us to reveal new and surprising features of many-body 
systems. In the Brout scheme, the mean ("molecular") field is spatially (z— ) independent. 
We find that there are no physical reasons for this assumption and extend the "mean-field 
concept." Within our approach, the so-called mean field is spatially dependent up to the 
moment when one should find the actual ground state. Then the fluctuation phenomena 
occur with respect to this ground state. In our approach the GL parameters and, hence, 
the ground state have a more precise evaluation. 

Our consideration has been performed for interactions J(|i — j|) of a quite general type, 
namely for all interactions that can be presented by the equality 



where Jq is an effective exchange constant, and z is an effective "coordination" number 
(number of interacting neighbours). In particular cases our results will be referred to the 
most common case of nn interactions (then n = 2D for simple lattices). For a simplicity, 
here we shall assume that IM is defined by a simple cubic (sc) lattice. The consideration 
can be easily expended to other types of regular lattices, as well as to irregular lattices 
with certain forms of quenched disorder, for example, random potential PP. We introduce 
the interaction radius by Rint ~ quantity, which is equal to the so-called zero- 

temperature correlation length (see Section 2). Let emphasize that our results can be 
rederived without difficulties in the presence of an external field {hi} conjugate to the 
lattice field {si}. 

In Section 2 we present a general approach to the treatment of fluctuation correlations 
at various length scales and discuss aspects of the usual theory that corresponds to the 
lowest order approximation of a perturbation expansion of cumulant type. In Section 3 we 
investigate the higher orders of the mentioned perturbation expansion and demonstrate 
several new features of the effective field theory. Our main results are summarized and 
discussed in Section 3.4-3.7. This report will be published in a more extended form jT3] . 

2. Present status of the effective field theory 

2.1. General scheme. The equilibrium free energy of IM as a function of the tempera- 
ture T and the field configuration h = {hi} is given by 




(2) 
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G{T) = -/5-Hn{Tre^^^(^)} , 



(3) 
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where = ksT, and the Trace is over the allowed lattice configurations s = {si}. 
Note, that the equilibrium values of the physical quantities are calculated as averages 
with respect to the statistical ensemble based on the Hamiltonian (P) and, in particular, 
the statistical averages of type {si...Sj) are obtained as derivatives of the partition sum 
whereas the irreducible averages {{si...Sj)) = ((sj — {si))...{sj — {sj))) are obtained as 
derivatives of the Gibbs free energy (jSl). 

Let us introduce the shift 

Si = + 6si , (4) 

where the lattice field 0j is an arbitrary (auxiliary) field configuration that is not neces- 
sarily associated with the averaged spin (sj) at site i, and the "fiuctuation" Ssi is merely 
the difference {si — (pi)- The identification of (pi with a statistical average (...) over the full 
Hamiltonian (0), or, with a statistical average (...)o corresponding to another ensemble as 
well as the interpretation of 6si as a fluctuation around (sj) (or (sj)o) may be a matter 
of further considerations. At this stage (pi and 6si are auxiliary variables which obey (jlj) 
and are not referred to concrete physical quantities. 

Following a standard procedure (see, e.g., Refs. [T1I2]) we obtain the following effective 
non-equilibrium free energy 

n{(P) = no{(P) + nfi(P), (5) 

with (p = {(pi}, 

^o(0) = IY1 ■^^i'^^'t'i - E 1^ [2ch(/3a,)] . (6) 

ij i 

Here 

ai = Y^ Jij(pj (7) 



is the "mean" (molecular) fleld, and 



nf{(P) = -p-Hn{^exp 
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is the "fluctuation" part. As usual, we shall often call the free energy (0) an "effective 
Hamiltonian." 

In (jHI), (...)o denotes a statistical average over an ensemble deflned by the auxiliary ("MF") 
Hamiltonian 

Ha{(p,s) = - ^QiSi ; (9) 

i 

the respective partition function and (nonequilibrium) free energy are given by Za{(p) = 
Tr [exp(— /37ia)], and Ga{(p) = — P~^\TiZa{(p) ■ For this simple ensemble, we have 

(si)o = th[/?a,(0)] . (10) 

The calculation of averages of type {si...Sj)o as well as "irreducible" averages of type 
{{Ssi...Ssj))o is also straightforward. These averages represent a form of fluctuation cor- 
relations, but they are just an auxiliary theoretical tool rather than real objects. In 
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contrast, the real objects, namely, (full) statistical averages (...) and ((...)) within the 
total Hamiltonian cannot be exactly calculated. 

2.2. Usual theory. In the framework of the usual theory, the "fluctuation" term Tif 
is ignored. This is the MFA. In the present general format of the theory, we have 
(self-consistency) equations of state {dT-C/d(t>i) = at fixed T (and h = {hi}) = 0) - one 
equation per a lattice vertex i: 

^J,,{0,-th[/3a,(0)]} = O, (11) 
j 

where (f) = {(pi}. It is easy to see that the number {z — 1) of nonzero terms {i 7^ j) in all 
sums (jlip is equal to the number of nonzero interaction constants acting on the site i: 
J{\i— j\) > for ao < \i — j\ < Rint] J{\i ~ j\ = for \i — j\ = R> Rint- The "equations 
of state" (|TT|) can be written in the simple form 

0, = th [Pai{^)] . (12) 

The equivalence of pij) and (fT^ can be easily proven for any number > 1. 
From (jlOj) and ()12|1 we obtain 

= th [Pai{^)] = 0, . (13) 

Thus in this quite general form of MFA (Hf ~ 0), we have: (sj)o = 0j, (...) = (...)o, 
and G(T,h) = ?^o(0) is the equilibrium free energy that corresponds to extrema (in- 
cluding minima) of ?i(0) ~ 'Hq{4>) - the non-equilibrium MF free energy given by the 
non-equilibrium (arbitrary) order parameter field 0j. By 0j from (jllj) we denote the equi- 
librium configuration of the latter; hereafter the "bar" of the equilibrium value (pi of (pi 
will be often omitted. Now one may perform a Landau expansion for small (pi in order to 
get other known forms of the MF theory. 

However, at the present stage of consideration we are interested in some more generality 
and for this reason we continue our discussion of the non-equilibrium free energy functional 

G(T/0) = niT/(P) ^ no{T/(P) (14) 

as given by Eq. (jUI). In (|14j). the dependence of the non-equilibrium free energy G on is 
denoted by "/0" because of the more special role of this variable, namely, the variation 
6G should be zero at thermal equilibrium and from this condition one obtains the possible 
thermal equilibria (p (alias "self-consistency condition"). 

The expansion of the log— term in ® up to fourth order in (pi yields the known result for 
the lattice 0^— theory of the Ising model: 

^0(0) = 2 ^ 'Jij(Pi4>j ~ 2 ^ JijJik(pj(pk + X] Ji'jJikJilJim(pj(pk(pl(pm ■ (15) 
ij ijk ijklm 

One may apply LWLA to this form of the theory but we shall follow a different path. 

2.3. Continuum limit. Using the rule 

^ /, = p y d'^xfix) j dxfix) , (16) 
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where p = (N/V) we can write Eq. (j2)) in the form 



J = pj d^RJ{R) 



(17) 



with R = \R\] R= (x — y). Note, that the factor p in (fTUj) - (PT|) can be avoided and this 
is the usual practice. In the latter case the physical dimension of the respective physical 
quantity is changed by a factor [V] ~ [L]^', for example, = [f{x)]. Of course, one 

may use both variants. 

In the continuum limit the effective Hamiltonian (jHI) corresponding to a zero external field 
[h = 0) takes the form 

T<o = ^j d^x4>{x)I[4>{x)] - pp-^ j d'^xln ch {(i'^I [<p{x)]] , (18) 

where 

I[cl>{x)]=p ! d''yJ{R)ct>{y) . (19) 



Now we apply LWLA in the form 

\<P{x)-<P{y)f <^\ct>{x)<P{y)\ (20) 
and under this assumption truncate the Taylor expansion 

to the second order m. R = {Ra}. 

Using the approximation (j2n|l . (fT9|l becomes 

I[^{x)] = J<P{x) + ^W^{x), (22) 
where J is given by (|T7jl . and 

J = P j d^RJ{R)R^ . (23) 

Now we substitute in (18), perform the expansion up to order 0^(cc) and to second 
order in V0(a;). Besides, we should keep in mind, that in expansion in powers of 0, 
we cannot distinguish between T and T^o except for the 0^— term where the difference 
between T and T^o should be kept only to the lowest nonvanishing order; in our case, this 
is the first order in (T — Tco): see, e.g., Ref. [I]. Following these notes, we perform at 
a certain stage of the calculation an integration by parts with the convenient boundary 
condition V0(a?) = and obtain the well known GL effective Hamiltonian 

7^o[0(^)]=P j d^xi^-^[Vct>{x)f + '^ct>\x)+u,ct>\x)^ , (24) 



6 



with 

co = ^J, ro(r) = A;s(T-r,o), «o = ^ • (25) 

Here T^o = (J/ks) and terms of order to = (T — Tcq) /Tcq -C 1 have been neglected in cq 
and uq 0, i.e. these two parameters are calculated at Tco [I1I2]- 

To clarify the result ()24j) - (j^^ we shall mention that J{R) for i? > is very small 
and can be ignored. Setting J{R) ~ Jo in ffTTj) . comparing the result with J = zJq from 
((2)), and noticing that p = 1/v ~ a^, one obtains z ~ {Rint/do)^ as should be, and 
J ~ JqR^^. In the same way one gets J ~ JR^nt- 

The results (E^ - (j2Sl) show that the energy of the spatially dependent configurations of 
the field depends on the interaction radius. The latter serves as a coherence (correlation) 
length of the field 0(a;); see the parameter cq given in ()25p. In order to clarify this 
point, let us consider the so-called zero-temperature correlation length [T], defined by 
eo = i{T = 0) = [-co/ro(0)]'/'. Using ^ and T^o = J/ke we obtain = RintV2D. 

2.4. Discussion. Note that the temperature range of validity of these considerations is 
to{T) <^ 1 and (kao) <^ n. The spatially dependent fluctuations correspond to a higher 
energy than the uniform configuration, and hence the latter contains the deapest (global) 
minima of the effective free energy nevertheless we have written N "equations of state" 
as a result of the minimization of the effective Hamiltonian. Now one can easily show that 
the variation of the Hamiltonian ()24|) with respect to the field will give again spatially 
dependent solutions but in the usual theory they are interpreted as "spatially dependent 
fluctuations" which, together with the uniform fluctuation 6(f) = (0 — (p) towards the 
stable state (f), are all fluctuations in the system in LWLA. But we know that this picture 
contains the approximation (sj) ~ (sj)o = (p- 

Our point of view is the following. The averages ()o should not be taken very seriously. 
They are an auxiliary tool in our consideration, and are not the final aim of our investiga- 
tion. We have used these averages only because they appear along our way of obtaining 
an effective Hamiltonian Ti. (or Tig in the lowest order of the theory) in which the sta- 
tistical degree of freedom 0j varies in a wide range of values (— oo < 0j < oo). This is 
the only consistent interpretation of our consideration performed so far. We can assume 
that up to now we have obtained nothing else but an effective free energy (effective lattice 
Hamiltonian) in terms of the new lattice field (pi. In a clear approximation, this effective 
model is given by (0); for the expansion in powers of (pi, see fITSj) . 

Another important aspect of our consideration is that the new Hamiltonian contains 
more ij— interactions than the original model (Q). The mathematical form of (jl5p gives 
a clear physical interpretation of these interactions: the first (pi(pj— term in the r.h.s. of 
()15|) describes an interaction that is quite similar to the original inter-spin interaction 
in whereas the second term of the same type in ()15|) describes an indirect two-site 
jfc— interaction that is mediated by the "z— spins". This means that the latter interaction 
JijJik has twice larger radius of action than the original Jij exchange. The four-point 
interaction given by the third term in the r.h.s. of ()15|) can also be described in the above 
style, introduced for the first time in this paper. 

Therefore, even at this early stage of consideration (7i ^ Tio) we see that the effective free 
energy exhibits effects of "statistical extension of the original inter-particle correlations 
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(interactions)". This "principle of growth of statistical correlations" is well known in the 
general phase transition theory. Here we show the concrete mechanism of the respective 
phenomenon and a systematic way of description of growth of statistical correlations. 
This point will become more clear from the results in the next Section. We shall see that 
the investigation along this path leads to a quite unexpected and intriguing picture. 

3. Beyond the standard theory 

3.1. Perturbation series. Here we consider the 0— contributions to the effective free 
energy (Hamiltonian) H(0) which are generated by the term 7ij(0). For obvious reasons, 
terms that are (/)j-independent will be omitted. 

It seems convenient to rearrange our theory by introducing the auxiliary variables Aj = 
(sjo — (pi), and cjj = (sj — Sio), where Sjo = (si)o- Then H can be written as an infinite 
perturbation series in powers of the (perturbation) Hamiltonian part 



ij ij 

The respective series can be presented in the form 

^ CO 

where 



2 

ij 1=1 



nf{<P)=^-^{S\{sA))oc, (28) 

(...)oc denotes the so-called connected averages P; for example, the average {Sj)Q{Sj)o 
is excluded from the connected {Sj)oc- For this cumulant (semi- invariant) expansion 
rules, similar to the Wick theorem in the perturbation theory of propagator type, are not 
available and one should perform the calculations with some caution. The term Ti^j^ is 
equal to zero, and this leads to a reduction of some infinite series in the next orders of 
the theory (/ > 1). 

3.2. Lowest order correction. Let us consider the first term in the r.h.s. of Eq. (j27j) 

and neglect all others. This yields 7i = (TYq + ^/) iii the form 

'H ~ -^^JijJik4'i4'k- — '^JijJikJji4'k4'i (29) 

ijk ijkl 

— ^ JijJikJilJim4>j4>k4>l4>m + — -hj J jkJilJimJin4>k4>l4>m4>n ■ 

ijklm ijklmn 

Performing this straightforward calculation one readily sees a very important property of 
the present theory, namely, that the first term in the r.h.s. of fll5|l is totally compensated 
by a respective counter term coming from the new contribution (~ AjAj) to the effective 
free energy. Besides, another new term twice compensates the second term in ()15p so that 
the term of type JJcpip now appears with a positive sign. The 0^— part of the effective free 
energy also undergoes a drastic change due to the A A— correction coming from Eq. (j28j) . 
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The same result can be obtained in a more general and, perhaps, more convenient way, 
if we add the AA— term in (P7j) to Tio from © before doing the expansion of Landau 
type. Then, within the same lowest order approximation for the series ()27|) we obtain a 
more general result for 7i, namely 

^ = -\Y. ^^jth(/3a,)th(/3a,)) + J^JUWa,) - Yl [2ch(/?ai)] . (30) 

ij ij i 

This form of Ti. clearly shows the lack of the simple Jij(f)i(j)j term describing the direct 
two-site exchange. In our further considerations we shall be faced only with inter-particle 
interactions (correlations), which are extended at distances larger than Rmt- In order to 
obtain one must expand the transcendental functions in (j^ : (/3aj) <C 1. 

The new forms ()29|) and (jHUj) of the effective free energy describes only indirect two-site 
interactions because the direct two-site interaction disappeared from our consideration. 
Now we are at a stage of description of correlations which extend up to 2Rij and larger 
distances. The tendency of the growing length scale of the interactions included in our 
effective free energy will be the main and unavoidable feature of our further consideration. 

Another very important feature of our findings is that in the simplest variant of the 
theory when the field is uniform (0j = 0) as well as in the GL variant in LWLA given 
by Eq. (j2D) the values of the Landau coefficients Cq, tq, and Uq do not change within the 
framework of accuracy of the effective 0^-field theory, where terms (Landau invariants) 
of order 0{t^) are ignored; t = {T — Tc)/Tc. This property seems to exist to any order of 
the expansion ()27|) in the limit of an infinite-range {z oo) initial interaction J^-. The 
1/z— corrections to the parameters (co,ro,uo) of the effective field theory are obtained 
from the (/ > 1)— terms in the series (j27|) . 

3.3. High-energy corrections of higher order. We have already mentioned that the 
(/ = l)-term in (j27p is zero. A calculation of the next two terms (/ = 2, 3) in (j27|) has 
been carried out in Ref. 0. Following Ref. [2j we can write Ti^p in the form 

^(2) _ P ^ Jfj \^ T T ^j^k N 

The result (jHH) gives the first (I/2;)— corrections to the parameters Co,ro, and Uq. Adding 
the result (j3ip to 7i(0) we obtain a form of the effective Hamiltonian 7i(0) which is more 
precise than the preceding ones. Let us write down this quite lengthy expression of the 
effective Hamiltonian: 

p2 p3 ^3 

Ti- = — t/jj Jjfc Jj70fc0« + — Jj^- Jjfc JiZ0fc0« ^ JijJikJjlJkm'Pl'Pm (32) 

ijkl ijkl ijklm 

+ '^vJikJilJim4>j4>k4>l4>m — (3'^ JinJij JikJilJnm4>j4'k4'l4>m 

ijklm injklm 

H ^ ^ ^ JinJipJnjJnkJnlJpm4^j4^k4^l4^m~\~ ^ ^ JinJipJijJikJnlJpm4^j4^k4^l4^m 

inpjklm inpjklm 

l3^ ^ — ^ 2 ^ — ^ 2 

— — JinJijJikJilJim4>j4>k4>l4>m — JinJijJikJilJim4'j4>k4'l4>m ■ 

injklm injklm 
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The result (|32j) shows that the indirect two-point interactions of type JJcfxp available in 
the effective Hamiltonian ()29|) do not exist in this higher accuracy of the theory. The 
interactions of type JJJ(j)(f) and JJJJcjxj) presented in the effective Hamiltonian (jH^ 
extend up to distances SRmt- The same is valid for the four point interactions included 
in (ISl- 

Within LWLA the Eq. (jS21) yields 



H=i5^(r + cA;2)|0(fc|2 + ^ E 0(fci)0(fc2)0(fc3)0(-fci - fc2 - fca) . (33) 
In (122)), 



k (^1,^2,^3) 



1 + ^ j Co, r= (^1 + ^ j fo, u= + (34) 
where fo = ksiT — Tc) is given by the "true" (renormalized) critical temperature 

Tc = T,o (l - ^) • (35) 

In deriving this lattice version of the effective Hamiltonian we have performed the lattice 
summations in (j32p in the reciprocal (fc)-space with the help of LWLA: J{k) ~ [J — Cok"^). 

The present results demonstrate a type of renormalization of the GL parameters (cq, tq, 
Uq) of the effective Hamiltonian due to l/2;-corrections. By a suitable choice of units, one 
of these parameters can be kept invariant, for example, equal to unity. Therefore, within 
a suitable normalization of the theory, the field (p^x) acquires a 1/z-correction as well |2]. 

The term Ti^p in ()27|1 has the form 

^(3) ^ .3 th{f3ai)th{f3aj) JijJgJji , . 

^ ~ ^ ch\l3ai)ch\l3aj) Q ^ c\i^{l3ai)di\(3aj)ch\l3ai) ^ ' 
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2'/ 



■ch^(/?a,)ch^(/5a,) 2 ^ ^ ch^(/3a,)ch^(/?a,) 

. . . A,-A;Anth(/ja,) 
3 5^ ch (/?a,) 

Let us consider the contribution of the term Ti^ to the quadratic (0^-) part of the 
effective Hamiltonian (jHHj) . In performing the calculations for both short-range (nn) and 
long-range {Rint ^ ^o) we have to evaluate again several lattice sums. Here we shall 
mention a particular sum, namely. 



^^JijJiiJji- (37) 



N 
iji 

which is equal to zero for nn interactions but gives a contribution (^ J'^ / z) for interaction 
radius Rmt > clq- Thus we introduce the following interpolation formula: E = k{z)/z, 
where < ^(2;) < 1 is an interpolation parameter which is supposed to be a smooth 
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V^{z) = 1 + — + — . (42) 



function of the coordination number z. The hmitine corresponds to nn 

interactions and the hmiting case k, = 1 corresponds to interactions of larger size. We 
suppose that the shape of the function k,{z) depends on details of the function J{R). 

Bearing in mind these notes we have calculated the quadratic contribution to the 

effective Hamiltonian which comes from the Hamiltonian part Tij given by (jHEI)- 
Here we present the results for the parameters Tco, c, and r which define the quadratic 
part H2 of H: 

n2{^) = lJ2^o'mm\\ (38) 
k 

with the (bare) correlation function 

G,\k) = V,{z) [T - T, + D2(^)T,p2fc2] _ ^39) 

Here 

T.= (l- — -— )r.„, (40) 

p,(,) = i + i±i^ + !i±^2^±15^!, (41) 

z 6z^ 

and 

2 + 3k , 14 + 30/t + 9K2 

The result ()35p for Tc{z) has been published for the first time |3j in case of nn interactions 
(k = 0); note, that there are errors in Ref. jSj for the functions T>i{z) and 1^2 (^)- 

The functions T>i{z) and 'D2{z) renormalize the field (f){k) and the vertices cq, and tq. 
For a total renormalization of the parameters of the theory up to the second order in the 
expansion we need to know the (I/2;) ^-correction to the vertex uq- We suppose that 
the calculation of this correction can be accomplished in the way described above; this 
"z-renormalization" has been discussed to first order in (l/z) in Ref. [2]. Here we wish 
to stress that within our extension of the theory the magnetic susceptibility Go(0) is Vi 
times smaller than the known MF susceptibility corresponding to Pi (00) = 1. 

The numerical coefficients in (j4(J|) - (j42j) indicate that real numbers z of nn like n = 
2,4,6 for simple lattices of spatial dimensionalities D = 1,2,3, respectively, give a good 
expansion parameter l/z. The 1/z-orrections are more substantial for the case of short- 
range interactions {Rmt ~ ^o), and one may suppose that for such interactions the (l/z)- 
series (HOJ - (jl^ are asymptotic; for the case of Tc, see a discussion of this tpoic in Ref. [1]. 
But even in case of asymptotic type of these series they may give more reliable results 
than the "bare" values (co,ro,Mo) of the Landau parameters; see arguments presented in 
Ref. 0. 

3.4. Critical temperature. The critical temperature given by pUj) can be compared 
with exact and reliable numerical (MC) results. Let us consider nn interactions (/t = 0). 
For one-dimensional (ID) IM, we know that = 0, MFA predicts = 2Jo//cs (in this 
case, z = 2D = 2), and (^01) yields = 5Jo/6kB- This is a quite good result for ID 
systems with very strong fluctuation effects. In 2D systems the fluctuations are not so 
strong and we find that reproduces the exact Onsager result (Tc = 2.27Jo//cb) with an 
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error of 22%, i.e., we have Tc = 35Jo/12A;b. For 3D systems our result is = 89Jo/18kB, 
whereas the best series analysis and MC results yield a difference of 9%: Tc = 4.5Jo/^b 
(see, e.g., Refs. fTHl ITT]). Our results seem quite rehable. 

3.5. Ground state. The l/^^-correction to the vertex uo has not been yet calculated 
although this calculation does not present difficulties. In this situation we shall give 
notion for the ground state energy by using the first order corrections to tq and uq- The 
equilibrium free energy per site / = [F/N) is given hy F = Ti and (jH!^ for the {k = 0)- 
Fourier amplitude 0(0), which minimizes /. Denote for convenience 0(0) = \/N(p ^ for 
the low-temperature ordered phase. From (33) we obtain / = — (r^/lGw) < whereas the 
"unrenormalized free energy is /o = (TCq/N) = — (rg/lGMo) < Thus, using (jMj) we have 
/(T) = (1 + 2/z)fo{T), which means that the new effective theory has a lower energy of 
the ordered phase. This is true also in the case of T = 0, where tq = — fc^^co = ~J- For 
the zero temperature (ground) state we have /o(0) = — 3J/4 and /(O) = —3(1 + 2/z)J/4. 
This result is also along the right direction because the MF theories (Tio) give unreliable 
high values of the ground state energy. The order parameter y^'^{T) = {—r/Au) is (1 — z) 
times smaller than the respective quantity v^o(^) — {~fo/^'^o) iii fhe usual theory based 
on 7^0- 

3.6. Effective interactions and phenomenon of growing of statistical corre- 
lations. We have explicitly shown the phenomenon of the growing length size of the 
interpaticle correlations in a classic system of interacting particles. To see this we have 
already introduced a new interpretation of the terms in the effective Hamiltonian (see 
Section 2.4). Let us consider the terms present in Ti. as terms describing concrete inter- 
site interactions. While the initial interaction Jij in IM ensures only two-site correlations 
(interactions), the effective Hamiltonians (fT3j) . (P^ . and (jH^ contain multi-site effective 
interactions. In contrast to the usual theory (fT3j) . where only extremely short-range effec- 
tive correlations are contained, the more precise effective Hamiltonians contain long-range 
two-site {4>i4>j) and four-site {(f)i4>j(f)k4'i) correlations, and all of these correlations are indi- 
rect, i.e. the correlation, for example, between two sites (ij) is mediated by one or more 
other sites {k, ...). A direct ( Jjj)-interaction is presented by the first term in the r.h.s. of 
p5|) but also the system exhibits two indirect correlations of type and 4'i(f)j4>k4'i given 
by the last two terms in the r.h.s. of (fT3j) . In the more precise variants of the theory, 
where a larger portion of the initial partition function has been calculated, the direct 
intersite interaction vanishes, and the particles are correlated only by indirect effective 
interactions. The length scale of these correlations grows in a monotonous way with the 
increase of the accuracy of the calculation, i.e. with the increase of the number I of the 
terms in the series fl27p . If we take the two-site correlations in the nn IM as an example, 
the maximal length of extension of these correlations in is ao, in ()15|) - 2ao, in (j^^ - 
3ao, in (j32j) - 4ao, i.e. [p — l)ao, where p is the maximal number of the summation indices 
in the terms of type in a given effective Hamiltonian. Surely, the number p tends to 
A^. This means that the most accurate effective field theory of many-body systems will 
correspond to (almost-) infinite range of correlations. 

The origin of these correlations is purely statistical. This effect is known and has both 
general formulation and application in many-body physics. Here we have established and 
described in details the concrete mechanism of this effect and, moreover, we have per- 
formed a demonstration of the remarkable picture of successive growth of the correlation 
length scale. 
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3.7. Final remarks. Obviously the (I/2;) -corrections are not the main point of discussion 
at the end of this paper of restricted length. Let us mention that the growth of the 
correlations discussed in Section 3.6 is not related with the l/^-corrections. It exists for 
any z, even in the "MF limiting case" of z —>■ 00, when the GL parameters keep their 
"initial" values (co,ro,Mo)- This effect follows from the fact that the terms in the initial 
Hamiltonian are compensated by the first "fluctuation" correction; see the first term in 
the r.h.s. of ()27p. At the next level of accuracy of the calculation, terms coming from 
the (/ = 2)-term in ()27p compensate the available terms and this process continues up to 
the incorporation of all particles in the correlation phenomenon; remember that the term 
corresponding to / = 1 is equal to zero. 

Here we emphasize that the terms in Tio - actually one of the most often used Hamiltonian, 
does not exist at all. They vanish just after the inclusion of the first correction to the 
usual theory; see AjAj-correction in ^I7\i . In place of these terms other terms with more 
complex structure come from the perturbation series ()27|). The outlined picture clearly 
indicates, that the terms which finally remain in the 0^-theory, are terms of type 

1 ^pM-ijM _ pMJM+l^^ M^N, (43) 

where obvious notations have been introduced; for example, jS^J"^ denotes the first term 
in the r.h.s. of (j29j) . The 0^ terms behave differently, because a lowest order term in 
J,namely, a term of type /5^J^</)^ appear at any step of development of the series (j27|) . 

At any stage of this surprising picture of the infinite series of successive modifications 
of the Hamiltonian both 0^- and 0^-terms keep their numerical coefficients equal to that 
in the usual GL Hamiltonian Tio- This is true within the whole scope of validity of the 
expansion in powers of 0. 

An important note, which should be emphasized is the following. While the sum 
is invariant with respect to the site i in regular lattices, the sum (jZj) depends on the 
site i. The reason is that the field configuration {0j} which takes part in ((7j) is not the 
equilibrium field. For the equilibrium field phi^^ the sum ([7j) will not dependent on the 
site i. This is consistent with the general notion that the equilibrium order in the volume 
of a homogeneous system in lack of effects of external fields, should be uniform. 

Our consideration justifies the GL fluctuation Hamiltonian. However, we have presented 
a new and quite surprising picture of the interparticle correlations, which reveals new 
remarkable properties of the GL theory. Apart from the l/2;-corrections to the GL pa- 
rameters, the structure of this theory is absolutely comprehensive as a tool for investiga- 
tion of large-scale correlation phenomena in many-body systems. We are certain that our 
findings have an application beyond the field of phase transitions. 
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